Multiscaling analysis of ferroelectric domain wall roughness 
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Using multiscaling analysis, we compare the characteristic roughening of ferroelectric domain 
walls in Pb(Zro.2Tio.s)03 thin films with numerical simulations of weakly pinned one-dimensional 
interfaces. Although at length scales up to Lma > 5 |^m the ferroelectric domain walls behave simi- 
larly to the numerical interfaces, showing a simple mono-affine scaling (with a well-defined roughness 
exponent £), we demonstrate more complex scaling at higher length scales, making the walls globally 
multi-amne (varying £ at different observation length scales). The dominant contributions to this 
multi-amne scaling appear to be very localized variations in the disorder potential, possibly related 
to dislocation defects present in the substrate. 

PACS numbers: 68.35. Ct, 77.80.Dj, 68.35.Dv, 68.37.-d 



Domain walls separating regions with different orien- 
tations of the order parameter in ferroic materials can be 
seen as elastic interfaces pinned by disorder, part of a gen- 
eral framework applicable to systems as diverse as growth 
surfaces and fractures pQ, flame and wetting fronts [2], 
and domain wall networks in the early universe [3]. In 
these systems, inter facial roughening with very different 
scaling properties [4] is predicted as a function of their 
(non)equilibrium state [5 j and the type of disorder [6H8]. 
From a practical viewpoint, understanding the behavior 
of the domain walls as elastic disordered systems allows a 
more accurate description of domain switching, growth, 
and stability, all key parameters for the successful im- 
plementation of (mult i) ferroic devices based on domains 
[9] [10] or domain walls [TTJ [12]. More broadly, ferroic 
epitaxial thin films provide an excellent model system for 
testing theoretical predictions, since parameters such as 
field, temperature, and defect density can be controlled 
over a wide range. 

Previous experimental roughening studies considered 
ferroic domain walls as equilibrated mono-affine inter- 
faces collectively pinned by weak, randomly distributed 
disorder [T3HT5] , for which particularly simple Gaus- 
sian scaling is expected. Whereas mono-affine interfaces 
present scaling properties (self-similarity) which can be 
described by a single scaling exponent at all length scales, 
in multi-affine systems, the scaling properties vary with 
the observation length scale, leading to a hierarchy of lo- 
cal scaling exponents [16 . Distinguishing between both 
cases through a multiscaling approach has proven useful 
in fracture studies, where multi-affine behavior was ini- 
tially reported [TTJ [18] , but subsequently attributed to 
finite size artifacts [19]. However, recent atomic force 
microscopy (AFM) measurements have shown that the 
disorder potential landscape in ferroelectric thin films 
can in fact be very complex, with strong individual pin- 
ning centers [20j [21] , and local variations in the disorder 
strength and universality class [22]. Potentially, these 
systems could therefore provide the sought-after experi- 



mental realization of more complex multi-affine scaling of 
interfacial roughening predicted by theory |8] , in which a 
rich and diverse disorder landscape - and possibly other 
interesting features - could be accessed. The multi-affine 
nature of such a system could be established by multi- 
scaling analysis [19] of individual domain walls or even 
different portions of a single domain wall, and directly 
compared with an ideal mono-affine model of weak col- 
lective pinning, in which the scaling behavior is exactly 
known. 

In this paper, we report on such a study of the rough- 
ening of ferroelectric domain walls in Pb(Zro.2Tio.s)03 
(PZT) thin films, compared to numerically simulated 
one-dimensional interfaces in random bond disorder 
(which preserves the local symmetry but not magnitude 
of the order parameter). Using multiscaling analysis, we 
demonstrate that the PZT domain walls can be seen as 
globally multi-affine interfaces composed of mono-affine 
segments separated by strong, highly localized variations 
in the disorder potential. The behavior of the mono- 
affine segments (with a roughness exponent £ = 0.57) 
is consistent with that of one-dimensional interfaces in 
random bond disorder up to a characteristic length scale 
Lma > 5 |Ltm, possibly related to the strong pinning or 
out-of-equilibrium effects of dislocation defects. 

Formally, for a roughened interface the geometrical 
fluctuations from an elastically optimal flat configuration 
at a given length scale r are quantitatively described by 
the probability distribution function (PDF) of relative 
displacements Au(r) — u(z) — u(z + r), whose character- 
istic scaling properties are reflected in the behavior of its 
central moments [4 : 

S„(r) = <|A«(r)|«>~r n S (1) 

where ( n are the associated scaling exponents for the nth 
moment, u(z) the transverse displacement along the lon- 
gitudinal coordinate z, (• • • ) an average over z, and ~ 
an average over disorder realizations when appropriate. 
One-dimensional interfaces at zero temperature equilib- 



2 



rium in uncorrelated disorder (weak collective pinning) 
are well described by a Gaussian PDF [7J [23j [24] , and 
are thus inherently mono-afflne. In this case, the rough- 
ness function B(r) = ^(Omono-affine ~ r<2< ^ suffices to 
fully characterize the scaling, with a single-valued expo- 
nent ("n = C^ n - Likewise, the displacement-displacement 
correlation functions, equivalent to the height-height cor- 
relation functions considered in fracture surfaces, 
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can be collapsed to a universal curve [19 . The specific 
value of the roughness exponent £ depends on the system 
dimensionality, the nature of the disorder, and the range 
of the elastic interactions. More complex disorder models 
[8], however, have shown multi-afflne behavior character- 
ized by an infinite set of individual scaling exponents ( n 
and a hierarchy of local roughness exponents. 

For the numerical study, simulations of a directed 
polymer model as a simple representation of a one- 
dimensional interface were performed on a discretized 
square lattice with an uncorrelated Gaussian noise distri- 
bution on each lattice site, using the solid-on-solid restric- 
tion \u(z + 1) — u(z)\ = ±1. The equilibrium zero tem- 
perature configuration was obtained using the transfer- 
matrix method with a droplet geometry, i.e. with one 
end pinned at the origin while the other end is free. 
The roughness exponent = 2/3 is exactly known for 
this model [25 , 26 j. To compare its statistical properties 
with those of the ferroelectric domain walls, we obtained 
10 4 equilibrium configurations from independent disor- 
der realizations for a system of size L = 2048, and deter- 
mined the Au(r) distributions for different length scales 
r. As shown in Fig. [T^l, these distributions are essen- 
tially Gaussian for intermediate length scales, with small 
deviations observed at very small and very large length 
scales. This observation is confirmed by the behavior 
of the displacement-displacement correlation functions, 
which collapse on a universal curve when renormalized 
by the Gaussian ratios = C^^/C^r), independent 
of r and ( [19] . Fig. [T]o shows such a collapse for orders 
n — 2-8, for the numerical interfaces in the r G [10 : 500] 
intermediate regime. We attribute the slight deviations 
once again apparent at very small and large length scales 
to finite size effects. Indeed, deviations at very small 
length scales are a known signature of the system mi- 
crostructure, in our case corresponding to the lattice dis- 
cretization. At very large length scales, as r approaches 
the system size, lack of statistics prevents sufficient aver- 
aging. 



From the numerical simulations, we extract £ a 



0.66 



from the B(r) averaged over the different disorder realiza- 
tions (Fig. [2^i), in excellent agreement with the theoreti- 
cal roughness exponent (j^ = 2/3. A comparable value 
of £ = 0.64 is obtained from the mean of the distribution 
of individual roughness exponents, shown in Fig. [2b. We 



FIG. 1. (a) PDF of the relative displacements for different 
length scales taken over 10 4 numerical disorder configura- 
tions in Gaussian units, compared with the Gaussian func- 
tion (solid line), a is the standard deviation, (b) Collapse of 
the Gaussian-normalized displacement-displacement correla- 
tion functions for orders 2-8. 



attribute the slight negative skewness of the histogram 
to the local solid-on-solid restrictions, constraining the 
roughness exponent to ( < 1, thus "compressing" the 
distribution to the right. In experimental systems pre- 
senting a symmetric £ distribution, the direct averag- 
ing method should provide a satisfactory estimate of the 
characteristic roughness exponent. 

An important, previously unremarked feature with sig- 
nificant consequences for the interpretation of both nu- 
merical and experimental results, is the large spread of 
the individual roughness exponent values. We emphasize 
here that mono-affine interfaces presenting a Gaussian 
PDF are expected to be fully described by a single- valued 
exponent. The large spread of the values extracted for 
C from individual numerical interfaces, with the distri- 
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FIG. 2. (a) Disorder- averaged roughness function with a 
roughness exponent Cavg = 0.66, and (b) roughness ex- 
ponents of all individual numerical interfaces, with mean 
value £ = 0.64, both in excellent agreement with theoretical 
£r§ = 2/3. (c) £avg and £ (shifted by 0.1 for visual clarity) 
for small numbers of interfaces, with the values for 10 4 config- 
urations indicated by long-dashed lines, and ±a by the short- 
dashed lines, (d) Converging FWHM of the £ distribution for 
small numbers of interfaces. 
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bution characterized by a full width at half maximum 
(FWHM) of nearly 30% of the mean value, therefore 
clearly shows that there exist inherent fluctuations of the 
scaling behavior from one interface to the other. To ob- 
tain a meaningful value of the roughness exponent, it is 
therefore crucial to average over a sufficiently large data 
set. From the rapid convergence of the FWHM of the 
distribution as a function of the number of interfaces in- 
cluded in the average (Fig.[2jl), we conclude that reason- 
able averaging can already be obtained with a few tens 
of different configurations. Moreover, the fluctuations of 
Cavg and ( for small numbers of interfaces are strikingly 
similar and equally rapidly bounded by a standard devi- 
ation a « 0.15 (Fig. [2^) . 

In parallel to these studies, experimental measure- 
ments of ferroelectric domain wall roughness were car- 
ried out on epitaxial ferroelectric PZT thin films, mon- 
odomain up-polarized as-grown, with the polarization 
axis perpendicular to the film plane [27]. Rectangular 
down-polarized domains with 180° domain walls were 
written using a positively biased scanning AFM tip, 
and imaged by piezoresponse force microscopy (PFM), 
as illustrated in Fig. [3^i, allowing the domain wall po- 
sition to be precisely determined. For each domain 
wall, we first computed the normalized displacement- 
displacement correlation functions C n (r)/R^ for orders 
n = 2-8. The resulting functions show roughnening up to 
a length scale r* = 100 nm (L* in [28]) and display behav- 
ior ranging from complete collapse to complete fanning, 
as shown for two domain walls less than 1 |Ltm apart in 
Fig. [3^-c. The observed saturation regime above r* may 
be attributed to the artificial writing process, limiting 
domain wall roughening at higher length scales [14] . Fan- 
ning between the correlation functions of different orders 
may of course be induced by instrumental resolution limi- 
tations [19] , and in our measurements the partial fanning 
observed for the smallest length scales can be attributed 
to exactly such size effects. However, the repeated ob- 
servation of radically different behaviors at intermedi- 
ate length scales between simultaneously-imaged domain 
walls emphasizes that these features are physically mean- 
ingful: the signature of the Gaussian and non-Gaussian 
nature of the underlying PDF, associated with mono- 
affine and multi-affine interfaces, respectively. 

We note, moreover, that on certain domain walls pre- 
senting fanning correlation functions when analyzed in 
their entirety, a collapse of the correlation functions could 
be recovered when only a segment of the domain wall 
was selected. In most cases, this selection corresponded 
to the exclusion of obvious local fluctuations, as demon- 
strated in Fig. [3] The exclusion of a visibly larger fluctu- 
ation in the middle of the right wall promotes the global 
non-Gaussian signature of Fig. ^ to the locally Gaussian 
signature of Fig. [3]i over the domain wall segment delim- 
ited by the short dashed box in Fig. [3^l, illustrating the 
strongly non-uniform and local character of the domain 



wall scaling properties. As expected for a multi-affine in- 
terface characterized by a varying roughness exponent ( 
at different length scales, distinct values of ( =0.78, 0.72 
and 0.66 were computed over segments of 1, 2 and 4 |Ltm 
centered on the middle fluctuation. 

Although the microscopic origin of such multiscaling 
behavior remains to be elucidated, and could be related 
either to the presence of strong pinning centers [8] , or to 
out-of-equilibrium behavior of the domain walls beyond 
a certain length scale [28 , our results convincingly show 
that the non- Gaussian nature of the relative displacement 
PDF, indicated by the observed fanning of the correla- 
tion functions, corresponds to a true multi-affine state. 
However, smaller segments of these domain walls remain 
mono-affine, presenting a Gaussian PDF. In other words, 
there exists a characteristic length scale below which do- 
main wall portions behave as mono-affine interfaces in 
weak collective pinning and above which they behave as 
multi-affine interfaces. We stress that, while roughen- 
ing is only observed up to r* in both cases, multi-affine 
behavior originates from highly localized individual fluc- 
tuations at least Lma apart along the length of the con- 
sidered domain wall segment. In our measurements, we 
see at most one of the local disorder fluctuations in a 
typical PFM domain wall image, giving a lower bound of 
Lma ~ 5 \xm for this length scale. One possible source for 
such localized fluctuations are dislocation defects prop- 
agating from the SrTiOs substrate, which can act as 
strong pinning sites. Transmission and scanning elec- 
tron microscopy studies of selectively etched single crys- 
tal SrTiOs revealed dislocation densities of ~10 8 cm -2 
[29] [30] , qualitatively in agreement with the observed dis- 
order fluctuations. Moreover, recent in-situ transmission 
electron microscopy observations of ferroelectric domain 
switching clearly show the pinning of domain walls by in- 
dividual dislocations [31 . The ability to non-invasively 
identify these very localized features at the domain walls 
could be especially interesting when coupled with an in- 
vestigation of their functional properties [32] . In particu- 
lar, given that oxygen vacancies [33, 34 increase domain 
wall conduction, and dislocation cores are associated with 
a very high presence of oxygen vacancies [35], exploring 
the link between mult i- affinity and domain wall current 
levels could be a promising research pathway. 

Finally, using multiscaling analysis, we were able to se- 
lect only mono-affine domain wall segments (on 43 out of 
63 domain walls imaged), for which a single- valued rough- 
ness exponent can be defined. For these domain walls, we 
constructed the roughness exponent distribution, shown 
in Fig. [4] with a mean value of 0.57 and a FWHM of 0.12, 
similar to the large distribution spread for the numerical 
interfaces, and identical to the C avg = 0.57, determined 
from B(r) (for 22 walls). 

This value is significantly higher than ( ~ 0.25 re- 
ported in previous studies of PZT domain wall rough- 
ness [14] . We believe a crucial difference between the two 




FIG. 3. 1.2 x 4.8 urn 2 PFM phase map showing two domain walls, with light and dark contrast corresponding to up and down 
polarization, respectively. For the full domain walls, the normalized correlation functions show either collapse (b) or fanning 
(c) between different orders. Collapse is recovered (d) when considering only the lower part of the right wall (dashed box). 
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FIG. 4. (a) Roughness exponents obtained on a set of 43 dif- 
ferent mono-affine domain walls with f = 0.57 =b 0.05. (b) 
Disorder-averaged roughness function for 22 domain walls, 
with similar roughness exponent Cavg = 0.57. 



experiments is the use of ultra-high vacuum, eliminating 
the screening effects of surface adsorbates. In both cases, 
the saturation of the roughness function B(r) clearly in- 
dicates that the artificial writing process limits domain 
wall roughening up to r*. At ambient conditions, r* is 
generally smaller or comparable to the film thickness. 
Coupled with independent studies of their dynamics, this 
yields an effective domain wall dimensionality of 2.5, in 
good agreement with theoretical predictions for elastic 
two-dimensional interfaces in random bond disorder with 
long range dipolar interactions [36]. However, in ultra- 
high vacuum r* extends to further length scales, allowing 
the interfaces to be considered as one-dimensional, with 
short-range elasticity. Moreover, progressively higher 
temperature thermal cycling of ambient-written domain 
walls with initial ( ~ 0.25 also promotes roughening, in- 
creasing ( values to ~ 0.5-0.6 [28 j. Although a Gaussian 
PDF cannot be taken to unambiguously imply that the 
corresponding domain walls are in fact in equilibrium, a 
scenario in which the mono-affine segments correspond 
to domain wall portions locally equilibrated with the un- 
derlying disorder landscape could agree with all the ex- 
perimental observations. 

The value of ( = 0.57 is close to previous reports for 
magnetic thin films (Cavg = 0-69 for 36 domain walls) [15] , 
for shallow periodic domains in ferroelectric single crys- 
tals (Cavg = 0.67 for 5 domain walls) [37 , and BiFeOs 



thin films (Cavg = 0.56 for 7 domain walls) [55] , all of 
which were taken to indicate one-dimensional domain 
walls weakly pinned by random bond disorder. However, 
none of these previous studies consider the possibility of 
deviations from mono-affine behavior, which could mod- 
ify the obtained ( values. 

In conclusion, our study of the scaling properties of 
numerical interfaces and ferroelectric domain walls has 
allowed us to directly compare ideal weak collective pin- 
ning with a real system presenting strong disorder fluc- 
tuations. Using multiscaling analysis as a tool to de- 
termine the (non) Gaussian nature of the PDF of relative 
displacements we find that the AFM- written ferroelectric 
domain walls show multi-afflne scaling, possibly related 
to the presence of dislocation defects. Between these local 
fluctuations smaller domain wall segments can be consid- 
ered as mono-affine, with an average roughness exponent 
C = 0.57. 
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